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Abstract 

We demonstrate that there appear finite-time future singularities in /(T) gravity with T being 
the torsion scalar. We reconstruct a model of /(T) gravity with realizing the finite-time future 
singularities. In addition, it is explicitly shown that a power-low type correction term (/3 > 1) 
such as a term can remove the finite-time future singularities in /(T) gravity. Moreover, we study 
/(T) models with realizing inflation in the early universe, the ACDM model. Little Rip cosmology 
and Pseudo-Rip cosmology. It is demonstrated that the disintegration of bound structures for 
Little Rip and Pseudo-Rip cosmologies occurs in the same way as in gravity with corresponding 
dark energy fluid. We also discuss that the time-dependent matter instability in the star collapse 
can occur in /(T) gravity. Furthermore, we explore thermodynamics in /(T) gravity and illustrate 
that the second law of thermodynamics can be satisfled around the flnite-time future singularities 
for the universe with the temperature inside the horizon being the same as that of the apparent 
horizon. 

PACS numbers: 04.50.Kd, 95.36.+X, 98.80.-k 
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I. INTRODUCTION 



A number of cosmological observations, e.g., Type la Supernovae cosmic microwave 



background (CMB) radiation {2, 3|, large scale structure (LSS) j^, baryon acoustic oscilla- 
tions (BAO) jsl, and weak lensing {g], have been implied the current accelerating expansion 
of the universe. Approaches to account for the late time cosmic acceleration are classified 
into two representative categories. The first is the introduction of unknown matters, Le., 
so-called "dark energy", in the framework of general relativity (for recent reviews, see J7^ 
The second is the modification of gravity such as f{R) gravity (for recent reviews, see |9K 

As a gravitational theory beyond general relativity, one could consider "teleparallelism" 
with the Weitzenbock connection, which has torsion T and not the curvature R defined by 



the Levi-Civita connection 



Q. 



In modern cosmology, in order to explain both inflation 



15| 



and the late time accelerated expansion of the universe, the teleparallel Lagrangian de nsity 



represented by the torsion scalar T has been extended to a function of T as /(T) (16|, |17 |. 
This idea is equivalent to the concept of f{R) gravity. In the recent literature, to check 
whether /(T) gravity can be an alternative gravitational theory to general relativity, its 



18 



the local Lorentz invariance 



22 



23|. 



various properties have been diversely explored 
non-trivial conformal frames and thermodynamics 

Moreover, it is known that if (phantom/quintessence) dark energy dominates the universe, 
in general there can appear finite-time future singularities, which have been classified into 
four types 2^. The finite-time future singularities in f{R) gravity 25[ have first been ob- 
served and those in various modified gravity Q, l^^] have also been investigated. Therefore, 
it is important to examine models of /(T) gravity in which finite-time future singularities 
can exist. 

In this paper, we concentrate on the two important theoretical features of f(T) gravity: 
the finite-time future singularities and thermodynamics in /(T) gravity. First, we explicitly 
reconstruct /(T) gravity in which the finite-time future singularities appear by following the 



procedure proposed in Refs. |28l-3o||. We also study a correction term to the models of /(T) 



gravity, so that such a term can remove the finite-time future singularities in analogy with 



f{R) gravity. Moreover, we explore the reconstruction of f{T) models which realize 



the 



examples of inflation in the early universe, the ACDM model. Little Rip cosmology 
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37| 



and Pseudo-Rip cosmology [38| . The Little Rip scenario is a kind of a mild phantom scenario 
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and considered in order to avoid a Big Rip singularity. On the other hand, the Pseudo- 
Rip model is an intermediate case between the cosmological constant and the Little Rip 
cosmology. In this model, the Hubble parameter asymptotically becomes constant as time 
goes to infinity, although for the Big Rip singularity, the Hubble parameter diverges at finite 
time, and in the Little Rip cosmology the Hubble parameter becomes infinity asymptotically 
as time goes to infinity. Furthermore, we examine whether the time-dependent matter 



instability in the star collapse 



occurs in /(T) gravity in analogy with f{R) gravity. 



This instability has recently been found in the framework of f{R) gravity, in addition to 



the well-known matter instability 



40|. Next, we explore thermodynamics in /(T) gravity. 



especially, near to the finite-time future singularities. In particular, we demonstrate that the 
second law of thermodynamics can be satisfied around the finite-time future singularities 
if the temperature of the universe inside the horizon is the same as that of the apparent 
horizon. We use units oi k-g = c = h = 1 and denote the gravitational constant SirG by 

= Svr/Mpi^ with the Planck mass of Mpi = G'^/^ = 1.2 x lO^^GeV. 

The paper is organized as follows. In Sec. II, we explain the fundamental formulations and 
basic equations in f{T) gravity. In Sec. HI, we investigate the finite-time future singularities, 
using analogy with f{R) gravity. In Sec. IV, we reconstruct /(T) gravity models with 
realizing the finite-time future singularities. We also study a correction term which can 
remove the finite-time future singularities. In addition, we reconstruct /(T) models with 
realizing inflation in the early universe, the ACDM model. Little Rip cosmology and Pseudo- 
Rip cosmology. The calculation of an inertial force which may lead to the dissolution of 
bound structures is done in Little Rip and Pseudo-Rip cosmologies for the Earth-Sun (ES) 
system. Furthermore, we discuss that the time-dependent matter instability in the star 
collapse can occur in /(T) gravity in analogy with f{R) gravity. In Sec. V, we explore 
thermodynamics in f{T) gravity. We demonstrate that the second law of thermodynamics 
can be satisfied around the finite-time future singularities. Finally, conclusions are given in 
Sec. VI. 
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II. /(T) GRAVITY 



A. Fundamental formulations 

Orthonormal tetrad components e^lx'^) are used in the teleparallehsm. An index A runs 
over 0, 1,2,3 for the tangent space at each point x'^ of the manifold. Their relation to the 
metric g'^'^ is described as g^^ = rjAB^^^u- Here, /i and u are coordinate indices on the 
manifold, which also run over 0,1,2,3, and e'^ forms the tangent vector of the manifold. 
The torsion T^,,„ and contorsion K^'^ „ tensors are defined as 

flu p 

T%.^e^A{d,et-d.e^) , (2.1) 
^ = -1 [T^'^ - T^'^p - ^'^) . (2.2) 

The teleparallel Lagrangian density is expressed by using the torsion scalar T, although in 
general relativity the Lagrangian density is described by the Ricci scalar R. The torsion 
scalar T is given by 

3^ = 5p"^n- (2-3) 

The modified teleparallel action describing /(T) gravity [l^ is as follows: 

^^^^ ' Cm] , (2.5) 



d x\e\ 

where |e| = det (e^) = \/—g and £m is the Lagrangian of matter. The variation of the 



action in Eq. (12. 5 p with respect to the vierbein vector field presents 16 1 



\d, {eS/n f - e^n,^, ''^f + S/'^d, (T) /" + -/aI = Ye^^T(^)/ , (2.6) 

where T^^^ is the energy-momentum tensor of all perfect fluids of ordinary matter, i.e., 
radiation and non-relativistic matter. 



B. Basic equations 

We take the four- dimensional fiat Friedmann-Lemaitre- Robertson- Walker (FLRW) space- 
time with the metric 

ds^ = Kpdx'^dx^ + f^dVL^ . (2.7) 



Here, f = a{t)r, = t and = r with the two-dimensional metric ha^ = diag(l, — a^(t)), 
a{t) is the scale factor, and dQ"^ is the metric of two-dimensional sphere with unit radius. 
In this background, we have g^^, = diag(l, — a^, — a^, — a^) and the tetrad components = 
{l,a,a,a). By using these relations, we find the exact value of torsion scalar T = — 6i7^ 
with H = a/a being the Hubble parameter, where the dot denotes the time derivative, d/dt. 

In the fiat FLRW background, the gravitational field equations can be written in the 
equivalent forms of those in general relativity: 

^' = y(pM + PDE) , (2.8) 
H = -YipM + Pm + PDE + Pde) , (2.9) 

where F = df/dT, F' = dF/dT, and pu and Pu are the energy density and pressure of 
all perfect fluids of generic matter, respectively. The perfect fluid satisfies the continuity 
equation pu + SH {pu + Pm) = 0. Moreover, the energy density and pressure of dark 
components can be represented by 

Pde = T^Ji , (2.10) 

Pbe = -:^{^J2 + Ji) , (2.11) 

with 

Ji = -T-f + 2TF, (2.12) 

J2 = {I- F -2TF')H . (2.13) 

Here, pde in Eq. fl2.1Up and Pqe in Eq. fl2.1ip satisfy the standard continuity equation 

Pde + 3i/ (pde + Pde) = . (2.14) 
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III. FINITE-TIME FUTURE SINGULARITIES IN /(T) GRAVITY 



A. Classification of the four types 



In the FLRW background (12. II) . the effective equation of state (EoS) for the universe is 
given by ^| 

2H 



Pcff 
-Peflf 



-Pcff 
Pes 



2H + 3H^ 



(3.1) 
(3.2) 
(3.3) 



Here, pcs and Pcfr correspond to the total energy density and pressure of the universe, 
respectively. When the energy density of dark energy becomes completely dominant over 
that of matter, one can consider wbe ~ Wcs- For H < {> 0), Wcs > —1 (< — 1), 
representing the non- phantom, i.e., quintessence (phantom) phase, whereas Wcs = — 1 for 
H = 0, corresponding to the cosmo logical constant. 



In Ref. 2Si, the finite-time future singularities has been classified into the following four 
types, (i) Type I ("Big Rip" [4l|): In the limit t — )■ tg, a — ^ oo, pcfr — ^ oo and |Pefr| — ^ oo. 
The case in which Ppff and Pcfr becomes finite values at t = tg 42| is also included, (ii) Type 



II ("sudden" [43|, |4J|): In the limit t tg, a as, Pca — ^ Ps and |Pcfr| — j- oo. (in) Type III: 
In the limit t ^ t^, a ^ a^, p^s — oo and |Pcfr| — oo. (iv) Type IV: In the limit t — ?• tg, 
a — )■ Og, pcfT — ^ 0, |Pcff| —J- 0, and higher derivatives of H diverge. The case in which 
and/or |Pcff| asymptotically approach finite values is also included. Here, tg, dsij^ 0) and ps 
are constants. 

It is important to mention that the Type I, i.e^"Big Rip" singularity, has recently been 



extended by Little Rip |3Ch37I| and Pseudo Rip 
the Type V ("u;") singularity, (v) Type V 
propagated (p.p.) curvature singularities 



scenarios. Furthermore, in addition to 



^"w" 



47|) singularity and earlier by parallel- 



48| have now been proposed. For the Type V 



("w") singularity, in the limit t — )■ tg, a — j- Og, pcS 0, |Pcff| — ?■ 0, and the EoS for the 
universe diverges. It should be cautioned that the Type V ("u;") singularity is similar to 
the Type IV singularity, namely, the Type V singularity does not lead to the divergence of 
physical quantities such as the scale factor, the energy density and pressure in the future. 
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but in the Type V higher derivatives of H do not diverge. This is the difference between 
the Type IV and Type V singularities. 

In what follows, since we examine the behavior of the universe around the finite-time 
future singularities, in which dark energy dominates over matter as Pde ^ Pm and Pqe ^ Pm 
in Eqs. ([23]) and (^M), we consider ~ Pde in Eq. P^s ~ ^de in Eq. f ETT]) . and 

WcS ~ W-DE = -Pde/PDE- 



B. Finite-time future singularities 



We consider t 
reference in Ref. 



le case in which the Hubble parameter H is represented by [the third 



28|] 



H 



hs 



{ts - ty ' 

hs 



for g > , 

for g < -1 , -1 < g < . 



(3.4) 

(3.5) 



{ts-ty 

Here, hs{> 0) and Hs{> 0) are positive constants, 0, —1) is a non-zero constant. More- 
over, tg is the time when the finite-time future singularity appears and only the period 
< t < ts is considered due to the fact that H should be real number. When t — )■ tg, for 
g > 0, both H hs {ts - t)"" and H ~ qhs (4 - 1)~^^^^^ become infinity. For -1 < g < 0, 
H is finite, but H becomes infinity. For g < — 1, but g is not any integer, both H and H are 
finite, but the higher derivatives of H can become infinity. It follows from Eq. (13.41) that 

hs 



a ~ as exp 



g-1 



hs 



{ts - t) 



for g 



-(9-1) 



forO<g<l, l<g, 



(3.6) 
(3.7) 



{ts-ty 

where Cs is a constant. 

It can be seen from Eq. fl3.6p that when t — tg? for g > 1, a — )■ oo, whereas for g < 
and < g < 1, a — Os. Moreover, it follows from Eqs. (13. 2p and (13. 4p that for g > 0, 
if — oo and therefore pcg = ?)H^ / — ?■ oo, whereas for q < 0, H asymptotically becomes 
finite and also p^s asymptotically approaches a finite constant value ps- On the other hand, 
from H ~ qhs{ts — t)"^'^'^^^ and Eq. ( 13. 3p we find that for g > —1, if — t- oo and hence 
Pcff = — (^2ii -|- 3ii^j /k^ — 7- oo. For g < —1, but g is not any integer, a, pcs and Pcs are 
finite because both H and H are finite, whereas the higher derivatives of H diverges. As a 
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TABLE I: Conditions for the finite-time future singularities to exist on q in the expressions of H 
in Eqs. (fOj) and ([33]) . pbe in Eq. ([230]) and Pde in Eq. (plT]) . and the behaviors of H and i7 in 
the Hmit of t — )■ tg- 

g(7^0, -1) H {t^Q H {t^t,) pDE i^DE 

g > 1 [Type I ("Big Rip") singularity] H ^ co H ^ oo 

or J2 / 

< g < 1 [Type III singularity] 
-1 < q <0 [Type II ("sudden") singularity] 
q < —1, but q is not any integer 
[Type IV singularity] 

derivatives 
of H diverge.) 



H ^oo 



H ^oo 
H ^oo 
H ^0 
(Higher 



J2 



result, the properties of the finite-time future singularities described by the expressions of 
H in Eqs. (13. 4p and (13. 5 P are summarized as follows: For q > 1, the Type I ("Big Rip") 
singularity, for < g < 1, the Type III singularity, and for — 1 < g < 0, the Type II 
("sudden") singularity. In addition, for g < — 1, but q is not any integer, the Type IV 
singularity appears. We present the conditions for the finite-time future singularities to 
exist on q in the expressions of H in Eqs. (13. 4p and (13. 5p . pde in Eq- (I2.10p and Pde in 
Eq. (12. lip , and the behaviors of H and H in the limit of t — )■ in Table HI 



IV. RECONSTRUCTION OF f{T) GRAVITY 

In this section, first we reconstruct /(T) gravity in which there appear the finite-time 
future singularities discussed in Sec. III^. Next, we examine a correction term removing the 
finite-time future singularities. 



^ In Appendix, we also describe a reconstruction method of f{T) gravity by way of using a scalar field 



through the extension of that of f{R) gravity 28I43C [ 
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A. Reconstruction 



By using Eqs. (12.101) and (12. lip , we find that the effective EoS for the universe at the 
dark energy dominated stage is written as 



DE 



4 (1 - F - 2TF') H + (-T - / + 2TF) 



PDE -T-f + 2TF 

As another description, from Eqs. f l2.1Qp and (12. lip we have 

Pbe = -PDE + HH, H) , 



(4.1) 



where 



2(1- F- 2TF') H 



(4.2) 



(4.3) 



Since T = —6H^, the form of /(T) is a function of H. It follows from Eqs. ( 13. 2p and (13. 3p 
that Peff = —PcS — 2H/k,'^. By comparing this equation with Eq. (14. 2p . we acquire the 
differential equation 

H+'^I{H,H)=0. (4.4) 
The substitution of Eq. (gj]) into Eq. (13]) yields 



H{F + 2TF') = . 
For H in Eq. (jMl), Eq. (gS]) reads F + 2TF' = because H ^ 0. 



(4.5) 



1. Power-law model 



As a form of f{T), first we take a power-law model, given by 



/(T) = AT^ 



(4.6) 



where 0) and 01(7^ 0) are non-zero constants. In this case, from Eq. (14. 5 p we have 



F + 2TF' = A (-6)""^ (2a - 1) //^(a-i) _ g . 



(4.7) 



In the limit of t — )■ ts, Eq. ( 14.71) has to be satisfied. From Eqs. ( 13. 4p and ( 13. 5p . we find that 
for g > (i.e., the Type I singularity [q > 1] and the Type III singularity [0 < q < 1]), 
a < 1, so that Eq. (14.70 can be satisfied asymptotically, whereas for g < (i.e., the Type II 
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singularity [—1 < g < 0] and the Type IV singularity [q < —1]), a = 1/2, in which Eq. (14. 7p 
is always satisfied. 

Furthermore, we state the meaning of the condition F + 2TF' = 0, which follows from 
Eq. fl4.5p . and another condition that the Friedmann equation fl2.8p can be satisfied, which 
may be interpreted as a consistency condition. Equations ( 12. 8 p and f l2.9p are expressed as 

T + Ji = , (4.8) 

J2-H = 0, (4.9) 

which are rewritten to 

-/ + 2TF = 0, (4.10) 

-F -2TF' = 0. (4.11) 

Here, we have used Eqs. fl212|) and fl27[3|) and if ^ for H in Eqs. ([33D and ([33]). It 
is clearly seen that Eq. (14.101) corresponds to a consistency condition and that the relation 
F + 2TF' = shown above is equivalent to Eq. (14. lip . Thus, the first condition is to satisfy 
the second gravitational equation (12. 9p . Moreover, from Eq. (I4.10p with Eq. (14. 6 p we find 
that the consistency condition becomes 

-/ + 2TF = A(-6)"(2a-l)if2" = 0. (4.12) 

In the limit of t — )■ ts, Eq. (I4.12p must be satisfied. For g > (i.e., the Type I singularity 
[q > 1] and the Type III singularity [0 < g < 1]), a < 0, so that Eq. (I4.12p can be satisfied 
asymptotically, whereas for g < (i.e., the Type II singularity [— 1 < g < 0] and the Type 
IV singularity [g < —1]), a = 1/2, in which Eq. (I4.12p is always satisfied. 
In addition, by using Eqs. (I2.12p and (I2.13P with (14. 6p . we obtain 

Ji = 6H^ [1-A (-6)"^' {2a - 1) ^ ^4 ^g) 

J2 = H[1-A (-6)"^' a {2a - 1) H^(^~^^] . (4.14) 

For the expression of H in Eqs. (13. 4p and (13. 5p . from Eqs. (I4.13P and (I4.14p we find that the 
conditions on Ji and J2 for the finite-time future singularities to exist in Table [T] is always 
satisfied. Thus, there can appear all the four types of the finite-time future singularities. 
We note that for a = 1/2, Eqs. KWf and fimp becomes = -T ^ and J2 = H 7^ 0, 
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respectively. Moreover, ii A = 1 and a = 1, this model in Eq. f l2.12p is equivalent to general 
relativity. 

It is very important to note that the conditions: [for g > (i.e., the Type I singularity 
[q > 1] and the Type III singularity [0 < g < 1]), a < 0, whereas for g < (i.e., the 
Type II singularity [—1 < g < 0] and the Type IV singularity [g < —1]), a = 1/2] derived 
in the above considerations are "necessary conditions" to produce the finite-time future 
singularities and not sufficient conditions. Indeed, if a < 0, the Type I singularity with 
g > 1 rather than the Type III singularity with < g < 1 appears because in the limit of 
t ts, both H and H with g > 1 diverge more rapidly than those with < g < 1. This 
originates from the absolute value of the power g, namely, the absolute value of g for the 
Type I singularity (g > 1) is larger than that for the Type III singularity (0 < g < 1). As a 
result, the Type I singularity is realized faster than the Type III singularity, and eventually 
the Type I singularity appears. Similarly, if a = 1/2, the Type IV singularity with g < — 1 
rather than the Type II singularity with —1 < g < occurs because in the limit of t — )■ ts, 
H ^ Hs and — ?■ wiht g < — 1 are realized more quickly than H ^ Hg and if — ?■ oo with 
— 1 < g < 0. This also comes from the absolute value of the power g, namely, the absolute 
value of g for the Type IV singularity (g < —1) is larger than that for the Type II singularity 
(—1 < g < 0). As a consequence, the Type IV singularity is produced faster than the Type 
II singularity, and accordingly the Type IV singularity appears. 

We also remark that in the Type V ("w;") singularity, a scale factor can be taken as {46 1 




n 



n - [2/ (3a)] 



n - [2/ (3a)] ^ ""^ 




(4.15) 



n — 1 

-1 y-' 

n-[2/{3a)]j 

where a and n are arbitrary constants. In the limit of t — > tg, H{t — )■ tg) — ^ and H{t — )■ 
tg) — !■ 0. On the other hand, the effective EoS for the universe = (1/3) (2gdcc — 1) — >■ 00. 
Here, gdcc = —da/b? is the deceleration parameter, which will again be redefined in Eq. fl4.3ip 
in Sec. IV C 2. Thus, in the limit of t — )■ ts-, since H{t — t- ts) = 0, from Eq. (14.51) we obtain 
F + 2TF' = 0. if we take a power-law model in Eq. (|46|) with A^O and a > 1, Eq. K7\i 
can be satisfied asymptotically because H{t — tg) = 0. As a result, if we take a power-law 

12 



model in Eq. (14. 6 p with A ^ and a > 1, the Type V ("w") singularity can appear. 



2. Exponential model 

Next, we examine an exponential model 

/(T)=Cexp(Ar) , (4.16) 

where 0) and 0) are non-zero constants. In this case, Eqs. (I4.10p and (14. lip 
becomes 

~ f + 2TF = C {-1 + 2XT) = 0, (4.17) 

-F - 2TF' = -ex (2AT + 1) exp (AT) = . (4.18) 

For the expression of H in Eqs. (13. 4p and (13.51) . in the limit of t — tg, both Eqs. (I4.17P and 
(I4.18P cannot be satisfied simultaneously. Thus, in an exponential model in Eq. (I4.16p there 
cannot apper the finite-time future singularities. 



3. Logarithmic model 

Next, we explore an logarithmic model 

f{T) = Dln{^T) , (4.19) 

where 0) is a non-zero constant and A(> 0) is a positive constant. In this case, 

Eqs. (I4.10p and (14. lip becomes 

-/ + 2TF = D[-ln(7T) + 2) = 0, (4.20) 
-F - 2TF' = ^ = . (4.21) 

For the expression of H in Eqs. (13. 4 p and (13. 5p . in the limit of t — )■ ts, both Eqs. (I4.20p and 
(I4.2ip cannot be satisfied simultaneously. Hence, in an logarithmic model in Eq. (I4.19p the 
finite-time future singularities cannot occur, similarly to the case of a exponential model 
in Eq. (I4.16P in Sec. IV A 2. Thus, in general the occurrence of the finite-time future 
singularities in /(T) gravity is realized in less cases than in f{R) gravity. 
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B. Correction term removing the finite-time future singularities 



It is known that in f(R) gravity, the addition of an term can cure the finite-time 
future singularities (see Ref. jol). Recently, it has also been demonstrated in Ref. j^^] that 
the addition of an term can remove the finite-time future singularities in non-local gravity. 
In this subsection, we investigate a correction term for the form of f(T) in Eq. (14. 6p so that 
the finite-time future singularities cannot appear. To execute this analysis, we explore an 
additional term of as function of T to the form of f(T) in Eq. f l4.6p so that for H in Eqs. (13. 4p 
and (13. 5p . the gravitational field equations (12. 8p and (12. 9p with Eqs. (I2.10p and (12. lip cannot 
be satisfied. 

As a straight forward procedure, we explore the case that the form of /(T) represented 
by Eq. (14. 6 P has a correction term fc{T), and analyze whether Eqs. (I4.10p and (14. lip can be 
satisfied or not. As an example, we choose a correction term fc{T) as 

f,{T) = BT^, (4.22) 

where 0) and 0) are non-zero constants. For /3 = 2, the correction term is similarly 
to that in f{R) gravity, i.e., a term. By combining Eqs. (14. 6 p and (I4.22p . the total form 
of /(T) including the correction term is expressed as 

/(T) = AT" + BT^ . (4.23) 

By substituting Eq. into Eqs. (I4:i0l) and (HlTll . we find 

- / + 2TF = A (2a - 1) + 5 (2/3 - 1) T'^ ^ , (4.24) 
-F - 2TF' = -Aa {2a - 1) T"~^ - 5/3 (2/3 - 1) T^"^ ^ . (4.25) 

From the considerations in Sec. IV A, we find that for g > (i.e., the Type I singularity 
[q > 1] and the Type III singularity [0 < g < 1]), /3 > 0, so that the second inequality in 
Eq. (I4.24P can be satisfied asymptotically, whereas for g < (i.e., the Type II singularity 
[— 1 < g < 0] and the Type IV singularity [g < —1]), /3 7^ 1/2, in which the second inequality 
in Eq. (I4.25P is always satisfied. As a result, if /3 > 1, for the Hubble parameter in Eqs. (13. 4p 
and (13. 5p . in the limit of t — )■ ts both of the gravitational field equations (12. 8p and (12. 9p 
cannot be satisfied. This means that a power-low type correction term with /3 > 1 can 
remove the finite-time future singularities in /(T) gravity. In Table UTt we describe necessary 
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TABLE II: Necessary conditions on the model parameters of a power-law model of f{T) in Eq. (j4.6p 
with realizing the finite-time future singularities, the emergence of the finite-time future singular- 
ities, and those of the correction term fc{T) = BT^ in Eq. (j4.'22p with removing the finite-time 
future singularities. 



g(/o, -1) 


Emergence 


/(T) = AT" 
(A / 0, a / 0) 


UT) = BTP 
{B^ 0,(3^0) 


q>l [Type I ("Big Rip") singularity] 


Yes 


a < 


P > 1 


< g < 1 [Type III singularity] 




Q < 


/3 > 1 


— 1 < g < [Type II ("sudden") singularity] 




Q = 1/2 


/3/1/2 


q < —1, but q is not any integer 


Yes 


a = 1/2 


/3/1/2 


[Type IV singularity] 









conditions on the model parameters of a power-law model of /(T) in Eq. (14. 6 p with realizing 
the finite-time future singularities, the emergence of the finite-time future singularities, and 
those of the correction term fc(T) = BT^ in Eq. (I4.22p with removing the finite-time future 
singularities. It is interesting to emphasize that a term, i.e., /3 = 2, which is the minimum 
integer to satisfy the condition /3 > 1, can remove all the four types of the finite-time future 
singularities in /(T) gravity. This consequence is the same as that in f{R) gravity. 

It is interesting note that for the case of the Type V ( "w" ) singularity, since in the limit of 
t — 7- ts H{t ts) = 0, if we choose a power-law type correction term /c(T) in Eq. ( ]4.22p with 
B 7^ and /3 < and substitute Eq. Km into Eqs. fHTOD and fl4TTll . we find Eqs. fOill 
and (14.251) . In other words, both of the gravitational field equations (12.81) and (12. 9p cannot 
be satisfied asymptotically. This means that a power- law type correction term /c(T) in 
Eq. (I4.22P with B ^ and /3 < can remove the Type V ( "w" ) singularity. 

C. Reconstructed models realizing cosmologies 

In this subsection, we explicitly present the reconstruction of /(T) models with realizing 
(a) inflation in the early universe, (b) the ACDM model, (c) Little Rip cosmology and (d) 
Pseudo-Rip cosmology. 
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1. Inflation in the early universe 



For generality, we consider power-law inflation. We suppose that the Hubble parameter 
is given by 

H=^, (4.26) 

where /iinf(> 1) is a constant larger than unity. It follows from Eq. fl4.26p that the scale 
factor is given by 

ait) = ainft'^'-f . (4.27) 

From Eq. f l4.27p . we flnd that d = Oinf/^inf (^inf — 1) t'^'"'"^ > 0, and hence power-law inflation 
occurs. In this case, since H = —h\ni/t^ 7^ 0, by using Eq. (14. 5 p we obtain the condition 
F + 2TF' = 0. 

As explained in Sec. IV A 1, for a power-law model /(T) = AT"' in Eq. (14. 6p . the 
conditions to be satisfled, which originate from the gravitational equations (12. 8p and (12. 9p . 
are given by Eqs. iK7} and (I4.12p . Therefore, if a < 0, in the limit of t 0, in Eq. (I4.26p 
diverges, so that Eqs. (14. 7p and (14.120 can approximately be satisfled in the very early 
universe. Moreover, if a = 1/2, Eqs. (14. 7p and (I4.12p can always be met. 

2. The ACDM model 

When we describe the ACDM model by the action in Eq. (12.50 . f(T) = T — 2A, where A > 
is positive cosmological constant, as is in general relativity. In this case, by substituting 
this form of f{T) into Eqs. (12.80 and (12. 9p . we have 

H' = ^, (4.28) 
H = 0. (4.29) 

Clearly, from Eqs. (^M) and fICTD . we flnd H = Ha = = constant, where we have 

deflned the Hubble parameter at the cosmological constant dominated stage as H\{> 0). 
Furthermore, the scale factor is expressed as 

a = OAexp {H\t) , (4.30) 

where aA(> 0) is a positive constant. 
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In the ACDM model, from Eq. (13. ip we find that the EoS is given by wbe = — 1 due to 
the fact that H is constant. Moreover, the deceleration parameter gdec, the jerk parameter 



j and the snark parameter s are defined by 



1 d-'^O, O-IN 

1 d^a 

^"3(,r-i/2)- (^•='3) 

By using Eq. fl4.30p . we obtain g^cc = — j = 1 and s = 0. 

The hmit on a constant EoS for dark energy in a fiat universe has been estimated as 
wde = —1.10 ± 0.14 (68% CL) [the second reference in Ref. [sl]. In addition, for a time- 



dependent EoS for dark energy with a linear form Wde(o) = w^deo + w;DEa (1 — o,) |50|, 



where wdeo and WDEa are the current value of wde and its derivative, respectively, the 
constraints have been analyzed as wdeo = —0.93 ih 0.13 and wdeu = — 0.41^o'7^ (68% CL) 
[the second reference in Ref. jsl]. Thus, the deviations of the values of (wde, Q'dec, s) from 
those for the ACDM model (—1, —1, 1, 0) show how the model is different form the ACDM 
model. In other words, we can use these four parameters as a observational test. 

We remark that if we consider the early universe, the model f{T) = T — 2A can lead to 
exponential inflation realizing de Sitter expansion of the universe, i.e., the Hubble parameter 
is given by 

H = i/inf = constant , (4.34) 
where Hi^f > 0. By using Eq. (14.341) . we have 

a{t) = ainf exp (iJinft) , (4.35) 

where ainf(> 0) is a positive constant. 



3. Little Rip cosmology 



Furthermore, we study Little Rip cosmology 30l ? -l36l]. which corresponds to a mild 
phantom scenario. The Little Rip scenario has been proposed to avoid the finite-time future 
singularities, in particular a Big Rip singularity. In this scenario, the energy density of 
dark energy increases in time with wde being less than —1 and then wde asymptotically 
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approaches wbe = — 1- However, such a scenario eventually leads to the dissolution of bound 
structures at some time in the future via the increase of an inertial force between objects. 
This process is called the "Little Rip" . 



As an example to realize Little Rip cosmology, we take the Hubble parameter as 



33| 



(4.36) 



where Hi^^{> 0) and ^(> 0) are positive constants. In this case, the scale factor a is 
expressed as 



LR 



exp (^t) 



a = flLR exp 

where aLR(> 0) is a positive constant. Moreover, from Eq. (13. ip we obtain 

2^ 



wbe = -1 



3H, 



LR 



■exp(-^t) . 



(4.37) 



(4.38) 



Since H = i^LR^ exp (^t) > 0, wde < —1, i-e., the universe is always in the phantom phase. 
In the limit of t — )■ oo, we find wbe —I and hence the Little Rip scenario can be realized. 

In the expression of wbe in Eq- (I4.38p . if we take ^ = Hq, at t = to ~ Hq^, we have 
wde = — 1 — ^Hq/ (SHi^Yi^). Here, to is the present time, Hq is the current value of the 
Hubble parameter given by Hq = 2.1h x 10~^^ GeV 5l| with h = 0.7 [the second reference 
in Ref. 3|, 52|], and e = 2.71828. By comparing this expression with the observational 
constraint on wbe = —1.10 ± 0.14 (68% CL) [the second reference in Ref. jsj], we find that 
if i^LR > [2Ho/{3e)] (1/0.24) = 1.50 x lO'^^ Q^y^ ^j^g current value of wbe in this Little 
Rip model is consistent with the observations. Here, we have used the fact that ^ = Ho and 
i^LR are positive values. 

By using Eqs. (l43T|) - fl4:33|) . ( K36\f and K37\f . we acquire 



i^LR exp (^t) 



J = l 



+ 3 



-f^LR [-f^LRexp (^t) J exp(^t) ' 
2ae + 3i/LRexp(et)] 



(4.39) 
(4.40) 
(4.41) 



For ^ = Hq, at t = to 



3i/LR [2e + 3i/LR exp (et)] exp (^t) ' 
Hq^, from Eqs. fl4.38p - p.4ip . we describe the expression of wbe 
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gdec, j and s at the present time to as 

WDE(t = to) = (4.42) 

Qdccit = to) = -1 - X , (4.43) 

j(t = to) = l + x(x + 3) , (4.44) 

s{t = to) = -^^^, (4.45) 

^ °^ 3 (2x + 3) ' ^ ^ 

with 

X = 7r^< 0.36, (4.46) 

where the second inequahty in Eq. fl4.46p follows from the observational constraint on Wj^e = 
-1.10±0.14 (68% CL) [the second reference in Ref. as x < (3/2) 0.24 = 0.36. As a result, 
if we take x ^ 1 enough for the deviation of the values of the four parameters (u^de, Qdccj, s) 
from those for the ACDM model (—1, —1, 1,0) to be very small, this Little Rip model can 
be compatible with the ACDM model. 

It follows from Eq. fl4.36p that in the limit of t — > oo, if diverges. For a power-law model 
/(T) = AT" in Eq. (gj]), if a < 0, in the hmit of t ^ cx), Eqs. g2D and KWf can be 
satisfied asymptotically. This is just an opposite case in the model with realizing inflation in 
Sec. IV C 1 because the limit in terms of t is the opposite direction. In addition, if a = 1/2, 
Eqs. (14.71) and (I4.12p can always be met, similarly to that in the case of inflation in Sec. IV 
C 1. 



4- Pseudo-Rip cosmology 



We also investigate Pseudo-Rip cosmology [3J, |38|. The above four cosmological models 
can be classified by using the behavior of the Hubble parameter as follows 38| . (a) power-law 
inflation: 

H{t) ^ oo , t^O. (4.47) 
(b) the ACDM model or exponential inflation: 



Hit) = H{to) . 



(4.48) 



(c) Little Rip cosmology: 



H{t) -> oo , t ^ oo . 



(4.49) 
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(d) Pseudo-Rip cosmology, which is also phantom asymptotically de Sitter universe: 



H{t) ^ Hoo<oo, oo. (4.50) 

Here, we consider t > to- Moreover, ifoo(> 0) is a positive constant. We also note that for a 
Big Rip singularity, H{t) — oo, t — > ts, as shown in Table [H As an example of a Pseudo-Rip 
model, we take 

H{t) = HpRtanhl-] , (4.51) 



.^0 

where i^pR(> 0) is a positive constant. In this case, the scale factor a is expressed as 

a = apRCOsh( — ] , (4.52) 



where apR(> 0) is a positive constant. From Eq. f l4.5ip . we find that H(t) is monotonically 
increasing function of t and H{t) — )■ iJpR < oo, t — ?■ oo. Thus, a behavior of H in the 
Pseudo-Rip cosmology in fl4.50p is realized. We also have H{t) = H^^J [tocosh^ (V'^o)] ^ 
0, t — 7- oo. This means from Eq. (12.91) that P — )■ —p in the limit of t — )■ oo. For a power-law 
model /(T) = AT" with a = 1/2 in Eq. (jM]), i.e., /(T) = A^T, Eqs. ([21]) and (Q can 
always be satisfied including in the limit of t — ?■ oo. 

For H in Eq. (I4.5ip . from Eq. (13. ip we find that the EoS is given by 

oto-npRsmh (t/to) 

From Eq. (I4.53p . we see that wd^ < —1, namely, the universe is always in the phantom 
phase, because H{t) = Hp^./ [to cosh^ (V^o)] > 0. In the limit of t — oo, we find wbe ^ — 1, 
similarly to that in Little Rip cosmology. 

It follows from wbe in Eq- (I4.53P that at t = to ~ Hq^, wbe = — 1 — 



v2 



[2Hq/ (SHpp)] 4/ (e — e ^)^ . In comparison with the observational constraint on wde 
-1.10 ± 0.14 (68% CL) [the second reference in Ref. Q], we find that if Hpn > 
{2Ho/3) 4/ (e - e'^f (1/0.24) = 2.96 x lO'^^ Q^y^ ^^e current value of wde in this Pseudo- 
Rip model is compatible with the observations. Here, we have used the fact that to ^ Hq^ 
and ifpR are positive values. 
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Moreover, by using Eqs. flOTjl - fOSj) . fOT]) and f02|) . we have 



(togpR)^tanh^ (t/to)-l 

fee = — 1 H n ^ , (4.54) 

(to//pR)'tanh2(t/to) 



l-(toiJpR.)'tanh' (t/to) 
(toi/pR)'tanh2 it /to) 
2 (toi/pR)'tanh2(t/to 



J = 1+ ^7TT X,:7:°^ (4.55) 



(4.56) 



2,toHp^ {toHpnf ia^-nh^ (t/to) + 2 ' 

At t = to ~ -f^cT^' from Eqs. f l4.53p - fl4.56p . we describe the expression of wde, Qdcc, j and 
s at the present time to as 

25 

WDE{t = to) = -1 - . > (4-57) 
3 smh 1 

5^ tanh^ 1-1 , , 

gdcc(t = to) = -1 + ^,^^^^2^ ' (4-58) 

. , , ^ 1 — 6^ tanh^ 1 , , . 

2 53 tanh^ 1-1 , , 

with 

5 = ^ < 0.497196, (4.61) 

-npR 

where the second inequahty in Eq. f l4.6ip follows from the observational constraint on wbe = 
-1.10 ±0.14 (68% CL) [the second reference in Ref. fl] as 5 < (3/2) 0.24 sinh^ 1 = 0.497196. 
Here, we use sinh^ 1 = 1.3811 and tanh^ 1 = 0.580026. As a consequence, we can take 
an appropriate value of 6 in order for the deviation of the values of the four parameters 
(wde, Qdccj, s) from those for the ACDM model (—1, —1, 1, 0) to be very small, so that this 
Pseudo-Rip model can be consistent with the ACDM model, similarly to that in the Little 
Rip model discussed in Sec. IV C 3. In Table IIIIl we display forms of H and /(T) with 
realizing (a) inflation in the early universe, (b) the ACDM model, (c) Little Rip cosmology 
and (d) Pseudo-Rip cosmology. 

In the expanding universe, the relative acceleration between two points separated by a 
distance / is given by Id/ a, where a is the scale factor. Suppose that there exists a particle 
with mass m at each of the points, an observer at one of the masses would measure an 
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TABLE III: Forms of H and f{T) with realizing (a) inflation in the early universe, (b) the ACDM 
model, (c) Little Rip cosmology and (d) Pseudo-Rip cosmology. 



Cosmology H f{T) 

(a) Power-law inflation H = h^ni/t , f{T) = AT"^ , 

(In the limit of t 0) /imf(> 1) a < or a = 1/2 



(b) ACDM model H = = constant , /(T) = T - 2A , 
or exponential inflation A > A > 

(c) Little Rip cosmology H = i^LR exp (^t) , f{T) = AT"" , 
(In the limit of t oo) i^LR > and ^>0 a<0ora = l/2 

(d) Pseudo-Rip cosmology H = Hp^ tanh (t/to) , /(T) = Ay/T 



inertial force on the other mass. The inertial force Fi^crt on a mass m is given by 31 



Fincrt = ml^ = ml{H + H^^ (4.62) 

= -ml J (pDE(a) + 3PDE(a)) = ml^ (^2pj,M + , (4-63) 

where in deriving the first equality in Eq. (I4.63P we have used Eqs. ( 12. Sp and ( 12. 9p . We take 
the present value of a as Oq = a{t = t^) = 1. We also provide that the two particles are 
bound by a constant force F^,. When -Fincrt(> 0) is a positive force and the amplitude is 
larger than that of Fb, the two particles become unbound and hence the bound structure is 
dissociated. In Pseudo-Rip cosmology, -Flnert is asymptotically finite. 

For a Big Rip singularity realizing H in Eq. (13. 4p with q > 1, by using Eq. (I4.62p we find 



inert 



mlh^ 



cx) , t^ts. (4.64) 



Here, the reason of the divergence is that H and H diverge in the limit of t — )■ tg- Moreover, 
for a Little Rip model in Eq. f l4.36p . we see that 

-Fincrt = "^^-^^LR + -f^LRexp (^t)] exp (^t) — ^ oo , t oo . (4.65) 

Similarly, the reason of the divergence is that H and H diverge in the limit of t — > oo. 
Thus, a phenomenon of "Rip" is produced by the cosmic accelerated expansion at a Big Rip 
singularity and in Little Rip cosmology. 
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On the other hand, for a Pseudo-Rip model in Eq. fl4.5ip . we obtain 



-^inert 

where 



mlH] 



PR 



to cosh [t/to] 



+ i7pRtanh" ( ^ 

tn 



7PR 

inert ,00 



pPR —mlfJ^ 
-'^ inert ,00 — ""-"PR • 



< 00 , t 00 , (4.66) 



(4.67) 



Therefore, -Finert is asymptotically finite. Here, the reason why Finert becomes a finite value 
of i^™t,oo ^he limit of t cx) is that H Hp^ and if 0. 



It is necessary for i^nSt oo be larger than the bound force F^^ = GMf^MQ/r"^ 



®-0 



4.37 X 10 GeV of the ES system in order that the ES system will be disintegrated and 



hence the Pseudo-Rip scenario can be realized. Here, 
Mr^ = 1.116 X lO^^GeV 



3.357 X 10^1 GeV 



7.5812 X 10^6 GeV"^ 



5l| and 



SI] are masses of Earth and Sun, respectively, and r, 



e-0 



lAU 



511 ] is the distance between Earth and Sun, i.e., the Astronomical 



unit. As an example, if we take m as the Earth mass m = M^, / is the distance between 



Earth and Sun / 



in order for -Fi^^t oo > ^h^^ by using Eq. fl4.67p . we find -ffpR > 



1.31 X 10 GeV. If this condition is met, the disintegration of the ES 



ffi-0 



system can occur much before arriving at de Sitter universe, so that the Pseudo-Rip scenario 
can be realized. In addition, if this constraint is satisfied, the current value of wbe in this 
Pseudo-Rip model is also consistent with the observations because the constraint on iipR 
from the current value of wbe is much weaker as Hp^ > 2.96 x 10~^^ GeV. 



D. Time-dependent matter instability and star singularity in f{T) gravity 

In modified gravity, the important property of the viability is the existence of the gravita- 
tional bound objects (stars, planets). It is checked via matter instability 40|. For instance, 
without such a property the relativistic star might not be formed because a corresponding 
singularity appears (for example, it is known that a singularity may appear for stars [53| 
in f{R) gravity). In a gravitating system with a time dependent mass density such as as- 



tronomical massive objects, the instability in f{R) gravity has recently been studied ^ 
Furthermore, the generation mechanism of the time-dependent matter instability in the star 
collapse has also been investigated. It has been demonstrated that the time-dependent mat- 



ter instability develops and consequently the curvature singularity could appear 54 [. In 



this subsection, by examining the process for the curvature singularity to be realized 



54| 



m 
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analogy with f{R) gravity, we discuss whether the time-dependent matter instabihty in the 
star collapse occurs in f(T) gravity. 

We examine a small region inside the star. We can regard this system as homogeneous 
and isotropic and hence the space-time is locally described by the fiat FLRW metric fl2.7l) . 
Here, the Hubble parameter H is negative because we are exploring the star collapse, so 
that the space-time can be shrinking. Since the region is shrinking, the energy densities of 
the matters automatically increase. In case of cosmology, as shown in Sec. HI B, all the 
four types of the finite-time future singularities in f{T) gravity can appear. If H diverges 
in the limit of t — tst, where tgt is the time when the time when the curvature singularity 
in the star appears, T becomes infinity because T = —6H^. From Table [H we see that for 
the Type HI singularity, T diverges, although the scale factor a is finite. This phenomenon 
can be applied to the star collapse. The energy density and the pressure from the matter 
are finite and therefore these can be neglected near the singularity because a asymptotically 
becomes a finite value. In this case, the Hubble parameter H is expressed as 



H — (4 68) 



where hst{> 0) is a positive constant. For the Type HI singularity, we have < g < 1. In 
the limit oi t ^ tst, T becomes infinite because H diverges. This means that the naked 
curvature singularity appears in the finite future. We note that in the above investigations, 
we have supposed that the region is almost homogeneous and isotropic. When these settings 
are adequate also near the singularity, the singularity simultaneously occurs in all of the 
region. The naked singularities occur densely in the region, even though the homogeneity 
and isotropy are broken. Moreover, the density of matter grows as in the region farther 
from the surface of the star, namely, nearer to the center of it. Accordingly, first the naked 
curvature singularity can appear near the center of the star. When the singularity produces 
the attractive force, the shrinking of the star proceeds more and more, whereas when the 
repulsive force is generated, eventually the explosion may happen. However, when the 
explosion happens, the sign of H has to change from negative to positive. The realization 
of this phenomenon seems to be difficult. 

On the other hand, in the cosmological context, when the Hubble parameter H is given 
by Eq. (13.41) with < g < 1, the Type HI singularity can appear. By using the initial 
condition for H and H, we can determine the values of tg (and hg) in Eq. (13. 4p and t^t (and 
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hst) in Eq. fl4.68p . The absolute values of H and H inside the star would be larger than 
those in the expanding universe because we are investigating the collapsing star. Hence, 
we find ts > tst- In other words, the curvature singularity in the star appears before the 
cosmological singularity. As a consequence, the time-dependent matter instability in the 
star collapse can occur in /(T) gravity, similarly to that in f{R) gravity. 



V. THERMODYNAMICS AROUND THE FINITE-TIME FUTURE SINGULAR- 
ITIES 



In this section, we discuss thermodynamics in f{T) gravity in order to examine whether 
/(T) gravity is a viable gravitational theory. In particular, by following the procedure in 



Refs. 
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55[, we examine whether the second law of thermodynamics can be verified around 



the finite-time future singularities . The fundamental connection between gravitation and 



56| (for recent reviews, see. 



thermodynamics was implied by Black hole thermodynamics 

e.g., js?!). In general relativity, by using the proportionality of the entropy to the horizon 



area, the Einstein equation was derived from the Clausius relation in thermodynamics 



This consequence has been app lied to more general extended gravitational theories 59 



58|. 
60|. 



A. First law of thermodynamics 



In Refs. 
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55 



61| , it has been shown that if the standard continuity equation (I2.14p in 



terms of dark component is satisfied, an equilibrium description of thermodynamics can be 
obtained. In the fiat FLRW space-time, the radius ta of the apparent horizon is written by 
= 1/H. The dynamical apparent horizon is determined by the relation h'^^^dafdpf = 0. 
The time derivative of = 1/H leads to —dfA/f\ = HHdt. By plugging Eq. (12. 8 p with 
this equation, we find the relation [1/ {AnG)] dfA = r\H (pt + Pt) dt. Here, pt = Pde + Pm 
and Pt = Pde + -Pm are the total energy density and pressure of the universe, respectively. 
In general relativity, the Bekenstein-Hawking horizon entropy is expressed as S = A/ (4G) 
with A = 4TTr\ being the area of the apparent horizon 56|]. By using the horizon entropy 
and the above relation, we acquire 



27rrA 



—dS = Anf^H (pt + Pt) dt 



(5.1) 



25 



The Hawking temperature Th = \K,sg\/ (2vr) corresponds to the associated temperature of 



the apparent horizon. Here, the surface gravity Kgg is represented by 
1 



62| 



1 



doc 
1 - 



Ta (1 - ^Wt) pt 



2tiG 



ta (pt - 3Pt' 



(5.2) 

(5.3) 
(5.4) 



where h is the determinant of the metric h^p and Wt = Pt/ Pt is the EoS for the total of 
energy and matter in the universe. It follows from Eq. (15. 4p that for Wt < 1/3, we have 
Ksg < 0. Thus, by substituting Eq. fl5.3p into Th = \i^sg\/ (2vr), we obtain 



Ta 



27lfA V '^Hta 
By combining Eq. (15. ip with Eq. (15. 5p . we acquire 

TndS = AnflH {pt + Pt) dt - 2Txf\ {pt + Pt) d^A 



(5.5) 



(5.6) 



The Misner-Sharp energy 63| is defined a.s E = ta/ (26*) = V pt with V = A7rf\/3 being the 
volume inside the apparent horizon. The last equality implies that E is equivalent to the 
total intrinsic energy. By using this equation, we have 



dE = -inr^H {p, + Pt) dt + iTir^PtdrA 



(5.7) 



We also define the work density [e^: W = - {1/2) [T'^^^'^^Kp + T^^^'^'^^Kfi) = 
(1/2) (pt — Pt)- Here, T^^^°'^ and T^^^'>°'^ are the energy-momentum tensor of matter and 
that of dark components, respectively. By using Eq. (15. 7p and the work density W, the first 
law of equilibrium thermodynamics can be described. 



TndS = -dE + WdV 



(5.8) 



As a result, an equilibrium description of thermodynamics can be realized. We also remark 
that from Eqs. (ESD, (Ell) and (O), we have S = Sn^Hf^ (pt + Pt) = (67r/G) (f/T^^ 
— {2t\: /G) H / {3H^) > 0. This means that in the expanding universe, the horizon entropy 
S always increases as long as the null energy condition pt + Pt > 0, i.e., H < 0, is met. 
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B. Second law of thermodynamics 



Next, we investigate the second law of thermodynamics in the equihbrium description. 
The Gibbs equation in terms of all matter and energy fluid is represented by T^dSt = 
d (ptV) +PtdV = Vdpt + (pt + Pt) dV. The second law of thermodynamics can be expressed 
as dSsuni/dt = dS/dt + dSt/dt > with S'sum = S + St, where St is the entropy of total 
energy inside the horizon. By using V = 47Tr\/3, Eqs. (12.91) . (15.51) and the relation S = 

M^m\ {pt + Pt) = 



(5.9) 




As a result, from the second law of thermodynamics with Eq. (15. 9 p we find the condition 23 1 
Y = - (aHT + = 12H (2H^ + if j > 0. From the behaviors of H and H in the limit 
of t — ts described in Table [H it can be seen that in the expanding universe, where H > 0, 
for all the four types of the finite-time future singularities, the relation 2H^ + H > can 
always be realized. It is also interesting to note that this relation is satisfied even in the 
phantom phase {H > 0). Thus, the second law of thermodynamics around the finite-time 
future singularities is always satisfied. However, it should be cautioned that at the exact 
singularity t = ts such as a Big Rip singularity, in which the scale factor a{t) diverges as 
O'it = tf.) = oo, a naive classical picture of thermodynamics might break down. We also 
mention that this result can be verified only provided that the temperature of the universe 
inside the apparent horizon is equal to that of the horizon 651 ]. 



VI. CONCLUSIONS 



In the present paper, we have reconstructed /(T) gravity in which finite-time future 
singularities appear. Furthermore, it has been demonstrated that a (/3 > 1) correction 
term, e.g., /3 = 2, to a model of /(T) gravity where the finite-time future singularities occur, 
can remove the finite-time future singularities, similarly to that in f{R) gravity. We have 
also explored non-singular /(T) models in which inflation in the early universe, the ACDM 
model. Little Rip cosmology and Pseudo-Rip cosmology are realized. It has been shown 
that the dissolution of bound structures for Little Rip and Pseudo-Rip cosmologies happens 
in the same manner as in gravity with corresponding dark energy fluid. Moreover, we have 
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considered that the time-dependent matter instabihty in the star collapse can occur in f{T) 
gravity, similarly to that in f{R) gravity. In addition, we have investigated thermodynamics 
in f{T) gravity and shown that the second law of thermodynamics can be satisfied around 
the finite-time future singularities, provided that the temperature of the universe inside the 
horizon is equal to that of the apparent horizon. 

It would be interesting to develop more complicated versions of f(T) gravity and study 
their cosmological applications. For instance, one can consider non-minimal coupling of 
/i(T) with electrodynamics, where /i(T) is a function of T, in analogy with that in non- 
minimal fi{R) gravity 66|. This may lead to the emergence of a domain- wall solution and 
variation of fine structure constant in non-minimal /i(T) gravity. From another side, it 
seems to be very interesting to generalize an /(T) theory in a consistent way so that one 
can include the presence of curvature in the Lagrangian. 

It should be emphasized that the illustration of the existence of the finite-time future 
singularities in f(T) gravity and the possibility of those removing due to an additional 
power- low term are nontrivial and significant. These consequences are also found in other 
alternative gravitational theories such as f{R) gravity. It is considered that the removal 
possibility of the finite-time future singularities can be one of the tests of a successful al- 
ternative gravitational theory to general relativity. Therefore, it is strongly expected that 
through these analyses of meaningful theoretical properties of modified gravitational theo- 
ries, we can obtain a successful alternative gravitational theory to general relativity which 
explains the cosmic accelerated expansion of the universe in a geometrical way. 
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Appendix A: Reconstruction method 



In this appendix, we explain the procedure of a reconstruction method of f{T) gravity 



by applying that of f{R) gravity 28|-|30| to /(T) gravity. The action of /(T) gravity with 
matter is given by Eq. fl2.5p . We can rewrite the action in Eq. (I2.5p with proper functions 
P{4>) and (5(0) of a scalar field (p to 

S = j d'x^^ iP{<l))T + g(0)) + j d'xCu {9,u, *m) . (Al) 

Since the scalar field does not have kinetic term, we may regard (p as an auxiliary scalar 
field. From Eq. (lAip . we derive the equation of motion of as = {dP{(f)) /d(f))T+ dQ {</))/ d(j). 
The substitution of = 0(T) into the action in Eq. (lAip leads to the expression of f(T) as 
/(T) = P{(f){T))T + Q{(f){T)). It follows from Eq. ^Ml that the gravitational field equation 
is described by 

rfP(0) d(j) 



e 



d, [eSD P{<P) - e\T^^,S/^P{<P) + S/^d, (T) ^— ^ 



+ \e-X + g(0)) = ye^T(^)/ . (A2) 

In the fiat FLRW background space-time with the metric in Eq. (12.71) . the components 
of (/i, z/) = (0,0) and (/i, i^) = («,j) = I,-- - ,3) in Eq. ( \A2h yield two independent 
differential equations in terms of P{(j){t)) and Q{(f){t)). In principle, by eliminating (5(0) 
from these two equations, we obtain an equation which P{(f){t)) obeys. We may take the 
scalar field as = t by redefining it properly. We express a{t) as a(t) = dexp {g(t)), 
where a is a constant and g{t) is a proper function of t. By using H = dg{(j))/ {d(f)), we 
represent the equation in terms of P{(j){t)) derived above to be a form described by P{(f){t)), 
the derivatives of P{(j)(t)) in terms of 0, dg{(f))/ {dip), and the derivatives of dg{(j))/ {dip) in 
terms of 0. Furthermore, from another equation we can also express Q{ip{t)) with P{ip{t)), 
the derivatives of P{ip{t)) in terms of 0, dg{ip)/ {dip), and the derivatives of dg{ip)/ {dip) in 
terms of 0. We derive the solutions of P{ip{t)) and Q{ip{t)) and substitute these solutions 
to f{T) = P{(p{T))T + Q{ip{T)), we acquire the expression of f{T) as a function of only T. 

Finally, we remark the following point. By redefining the auxiliary scalar field as 
= ^{(p) with a proper function $ and defining P{(p) = P{^{ip)) and Q{(p) = Q{^{ip)), the 
new action S = J d^x^f{T)/ {2k'^) + / d^xCu (^^.^, ^m), where f{T) = P{^)T + Q{^), 
is equivalent to the action in Eq. (lAip . This is because f{T) = f{T). Since ip is the inverse 
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function of $, by using = 4>{T), ip can be solved with respect to T as <y9 = <^{T) = 
$^^(0(T)). Accordingly, there exist the choices in as a gauge symmetry and therefore 
can be identified with time t as = t. We can interpret this fact as a gauge condition 
which is equivalent to the reparameterization of = (f){ip) 67|]. As a result, we can find 
the form of f{T) by obtaining the relation t = t(T). In addition, regarding the relation 
between H and T in f(T) gravity and that between H and R in f{R) gravity in the flat 
FLRW background space-time, we note that for /(T) gravity and f{R) gravity, we have 
T = —6H^ and R = 6 (^H^ + H^, respectively, and therefore in comparison with f{R) 
gravity, in /(T) gravity the torsion scalar T depends on only H, although in f{R) gravity 
the scalar curvature R depends on both H and H. 
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